HAW 

HAMBURG 


Projekt 


Überprüfung einer einfachen Kopfrechenmethode 
zur Umrechnung der Fluggeschwindigkeit 
von CAS in TAS 


Autor: Dennis Lucht 


Betreuer: Prof. Dr.-Ing. Dieter Scholz, MSME 

Abgegeben: 30.04.2019 


Fakultät Technik und Informatik 
Department Fahrzeugtechnik und Flugzeugbau 





1 


DOI: 

https://doi.org/10.15488/941 1 

URN: 

https://nbn-resolving.Org/um:nbn:de:gbv: 18302-aero2019-04-30.018 

Associated URLs: 

https://nbn-resolving.Org/html/urn:nbn:de:gbv: 18302-aero2019-04-30.018 

© This work is protected by Copyright 

The work is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 

International License: CC BY-NC-SA 

https://creativecommons.Org/licenses/by-nc-sa/4.0 



Any further request may be directed to: 

Prof. Dr.-Ing. Dieter Scholz, MSME 
E-Mail see: http://www.ProlScholz.de 

This work is part of: 

Digital Library - Projects & Theses - Prof. Dr. Scholz 
http://library.ProfScholz.de 

Published by 

Aircraft Design and Systems Group (AERO) 

Department of Automotive and Aeronautical Engineering 
Hamburg University of Applied Science 

This report is deposited and archived: 

• Deutsche Nationalbiliothek (https://www.dnb.de) 

• Repositorium der Leibniz Universität Hannover ( https://www.repo.uni-hannover.de) 

• Internet Archive (https://archive.org) 

Item: https://archive.org/ details/T extLucht.pdf 


This report has associated published data in Harvard Dataverse: 
https://doi.org/ 10.7910/DVN/VU7TB8 















2 


Kurzreferat 

Zweck - Von Piloten werden u.a. die sogenannten "Manual Flying Skills" gefordert. Dabei 
muss der Pilot in der Lage sein (ohne Autopiloten) nach grundlegenden Instrumenten zu flie¬ 
gen. Dafür sind neben Geschick auch Faustformeln erforderlich. Die Faustformeln müssen da¬ 
für verlässlich sein. Der Inhalt dieser Arbeit beschäftigt sich exemplarisch mit einer Faust¬ 
formel zur Umrechnung der kalibrierten Fluggeschwindigkeit (Calibrated Airspeed, CAS) in 
die wahre Fluggeschwindigkeit (True Airspeed, TAS). 

Methodik - In Excel und Matlab werden die Ergebnisse aus den Berechnungen der Faust¬ 
formel mit dem Ergebnis einer exakten Berechnungsweise anhand flugmechanischer Fonnein 
verglichen. Dabei wird die Flughöhe und Fluggeschwindigkeit variiert. Es werden die Abwei¬ 
chungen ermittelt und in Diagrammen zwei- und dreidimensional visualisiert. 

Ergebnisse - Die zu prüfende Faustformel liefert in dem für Sie vorgesehen Anwendungsbe¬ 
reich hinreichend genaue Ergebnisse mit Abweichungen unter 5 %. Dabei nehmen die Ab¬ 
weichung zu, umso weiter die Parameter (Höhe und Geschwindigkeit) von typischen Reise¬ 
flugbedingungen entfernt sind. 

Bedeutung in der Praxis - Piloten können bede nk enlos auf die in dieser Arbeit geprüfte 
Faustformel zurückgreifen und kommen so mit überschaubarem Kopfrechenaufwand auf rela¬ 
tiv genaue Ergebnisse. 

Wert - Diese Arbeit zeigt, wie mit mäßigem Zeitaufwand in Excel eine Faustformel über ei¬ 
nen gesamten Bereich geprüft werden kann. Das Vorgehen kann auf weitere Faustformeln 
übertragen werden, sodass sich ein Pilot sein "Kniebrett" mit verifizierten Faustformeln füllen 
kann. 



I =- HAW 
— HAMBURG 


DEPARTMENT FAHRZEUGTECHNIK UND FLUGZEUGBAU 

• • 

Überprüfung einer einfachen Kopfrechenmethode 
zur Umrechnung der Fluggeschwindigkeit 
von CAS in TAS 

Aufgabenstellung für ein Projekt 

Hintergrund 

Von Piloten werden u. a. die sogenannten "Manual Flying Skills" gefordert 
(https://www.skvbrary.aero/index.php/Pilot Handling Skills) . Dabei muss der Pilot in der 
Lage sein nach grundlegenden Instrumenten zu fliegen. Dafür sind neben Geschick auch 
Faustformeln erforderlich. Die Faustformeln müssen dafür verlässlich sein. In dieser Arbeit 
geht es exemplarisch darum, eine Faustformel zur Umrechnung der Kalibrierten Flugge¬ 
schwindigkeit (Calibrated Airspeed, CAS) in die Wahre Fluggeschwindigkeit (True Airspeed, 
TAS) zu überprüfen. Diese Faustformel wird von Piloten angewendet, die schnellere Flug¬ 
zeuge fliegen und wird in der Vorlesung "Flugmechanik" an der HAW Hamburg erklärt 
(http://www.fzt.haw-hamburg.de/pers/Scholz/materialFMl/l-19JPG) . In einem Beispiel 
(http://www.fzt.haw-hamburg.de/pers/Scholz/materialFMl/l-20JPG) wird gezeigt, dass die 
Berechnung (in diesem Beispiel) gut funktioniert. 

Aufgabe 

Die Aufgabe dieser Projektarbeit ist es, die Faustformel zur Umrechnung von CAS in TAS 
auf ihre Gültigkeit im gesamten Flugbereich zu überprüfen. Beginnen Sie bei den relevanten 
Gleichungen der Flugmechanik. Berechnen und visualisieren Sie den relativen Fehler abhän¬ 
gig von Flughöhe und Fluggeschwindigkeit. 

Die Ergebnisse sollen in einem Bericht dokumentiert werden. Es sind die DIN-Normen zur 
Erstellung technisch-wissenschaftlicher Berichte zu beachten. 
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Verzeichnis der Begriffe und Definitionen 

CAS 

"Die CAS (kalibrierte Fluggeschwindigkeit) ist die Geschwindigkeit, die direkt durch das 
Pitot Static-System des Flugzeugs gemessen wird. Sie entspricht der IAS, korrigiert um 
den Statikdruckquellenfehler. (...) In modernen Flugzeugen werden die Fehler durch einen 
Air Data Computer kompensiert, so dass die dem Piloten angezeigte Geschwindigkeit (IAS) 
der CAS entspricht." (Scheiderer 2008) 

Druckhöhe 

"Engl.: Pressure Altitude. Gibt für einen gemessenen Luftdruck an, welcher Höhe er in der 
Standardatmosphäre entspricht. Die Druckhöhe ist also die Flughöhe, die ein barometrischer 
Höhenmesser bei QNE-Einstellung anzeigt." (Klußmann 2007) 

EAS 

"Abk. für Equivalent Air Speed. Im Deutschen auch als äquivalente (Flug-) 
Geschwindigkeit oder Äquivalenzgeschwindigkeit bezeichnet. Die EAS entspricht der 
am Fahrtmesser im Flugzeug angezeigten Geschwindigkeit (IAS), die um den Einfluss 
der Kompressibilität der Luft bereinigt wurde. Diese Korrektur ist ab 
Fluggeschwindigkeiten nahe Mach 1 relevant." (Klußmann 2007) 

Flugfläche 

"Auch Flight Level genannt und mit FL abgekürzt. Bezeichnung für Flächen 
konstanter Druckhöhe, die als Höhenbänder der vertikalen Regelung des Flugverkehrs 
dienen. Flugflächen sind stets Vielfache von 500 ft und werden als Druckhöhen in Vielfachen 
von 100 ft angegeben. Die Angabe FL 65 steht also für eine Druckhöhe von 6500 
ft." (Klußmann 2007) 

Geopotenzielle Höhe 

"Die geopotentielle Höhe H ist eine Hilfsgröße, mit der die auf die Masse bezogene 
potenzielle Energie eines Fluidelements unter Berücksichtigung der Höhenveränderlichkeit 
der Erdbe-schleunigung beschrieben werden kann (...) Benutzen wir also anstelle der 
tatsächlichen Höhenkoordinate z die geopotentielle Höhe H, so können wir mit konstanten 
Normerdbeschleunigun go rechnen" (Kümmel 2007) 
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IAS 

"Abk. für Indicated Air Speed. (...) Bezeichnet die Fluggeschwindigkeit, die von einem 
Fahrtmesser angezeigt wird. Die Fluggeschwindigkeit wird vom Fahrtmesser indirekt über die 
Messung des Staudrucks ermittelt. (...)“ (Klußmann 2007) 

ISA 

"ISA bedeutet International Standard Atmosphere oder ICAO Standard Atmoshere. Eine 
Standard Atmosphäre ist eine hypothetische vertikale Verteilung der atmosphärischen Tempe¬ 
ratur, des Luftdrucks und der Luftdichte entsprechend einer nationalen oder internationalen 
Vereinbarung. Dies soll es beispielsweise ermöglichen, Höhenmesser zu kalibrieren, Flug¬ 
zeuge nach einheitlichen Kriterien zu konstruieren und Leistungskurven von Fluggerät zu 
entwickeln. “ (Schmidt 2000) 

Machzahl 

"Angabe der Fluggeschwindigkeit in Vielfachen der momentanen Schallgeschwindigkeit. Die 
momentane Schallgeschwindigkeit ändert sich dabei unter anderem in Abhängigkeit von der 
Dichte und der Temperatur der Luft. (...) Die Einheit der Machzahl ist das Mach.“ (Kluß¬ 
mann 2007) 

„Normatmosphäre ist ein Begriff aus der Luftfahrt und bezeichnet idealisierte Eigenschaften 
der Erdatmosphäre. Die internationale Standardatmosphäre (engl. International Standard At¬ 
mosphere, ISA) ist von der International Civil Aviation Organisation (ICAO) definiert wor¬ 
den. Sie stellt eine Atmosphäre dar, bei der die Größen Luftdruck, Lufttemperatur, Luftfeuch¬ 
tigkeit sowie Temperaturabnahme je 100 m Höhenstufe Werte haben, die ungefähr gleich den 
auf der Erde herrschenden Mittelwerten sind. “ (Nguewo 2018) 

TAS 

"Abk. für True Air Speed. Im Deutschen auch als wahre Eigengeschwindigkeit bezeichnet. 
(...) wirkliche Geschwindigkeit gegenüber der umgebenden Luft. (...) wobei zur Umrech¬ 
nung die Luftdichte in Meereshöhe verwendet wird. (...) Darüber hinaus berücksichtigt die 
TAS auch Einfluss der Kompressibilität der Luft, die ab Fluggeschwindigkeiten nahe Mach 1 
relevant ist. “ (Klußmann 2007) 

TAT 

"Die mittels Fluginstrument angezeigte Temperatur der umgebenden Luft, erhitzt durch die 
adiabatische Kompression (ram effect). TAT = SAT (1,2-M 2 )“ (Scheiderer 2008) 
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1 Einleitung 

1.1 Motivation 

Heute wird ein Großteil der Strecken in der zivilen Luftfahrt mittels Autopiloten geflogen. 
Dennoch sind in einigen Situationen die sogenannten "manual pilot skills" von besonderer 
Bedeutung. Damit ist der Blindflug ohne Autopiloten gemeint, was sich aber nicht nur auf das 
Geschick des Piloten beschränkt, sondern auch beinhaltet, dass der Pilot mit Rohdaten (von 
grundlegenden Instrumenten) fliegen kann. Dafür sind nebst diversen Finessen auch Faust¬ 
formeln erforderlich. Hierbei ist es von entscheidender Bedeutung, dass der Pilot sich auf die 
Ergebnisse der verwendeten Faustformeln verlassen kann. 

Ein Ergebnis der unzureichenden sogenannten "blind flying skills" bei Piloten ist die An z ahl 
der lost of control in flight (LOC-I) Unfälle. Fast alle dieser Geschehnisse ereigneten sich in 
der Nacht oder bei „Instrument Meteorological Conditions (IMG)“, wenn Piloten gezwungen 
waren, die Kontrolle zu übernehmen. 

Von der FAA wird angenommen, dass der kontinuierliche Einsatz von Autoflight-Systemen 
in dem zunehmenden Unvermögen resultiert, das Flugzeug schnell aus einem unerwünschten 
Zustand heraus zu stabilisieren. Ergo sollten Piloten zusätzlich und regelmäßig entsprechende 
Kenntnisse aufarbeiten, um das Risiko zu minimieren. Hierbei ergibt sich ein Kosten-Nutzen- 
Problem, denn das statistische Risiko tödlicher Unfälle ist aufgrund des Erfolgs moderner 
Cockpit-Technologie sowie einer präzisen und zuverlässigen Flugdeckautomation gering, so 
dass Piloten kaum manuell mit dem Flugzeug umgehen müssen. Aber wenn sie die Kontrolle 
übernehmen müssen, zeigen Statistiken, dass nicht garantiert ist, dass Piloten gute Fertigkei¬ 
ten im Hinblick auf die manuelle Flugsteuerung vorweisen können. Je weniger sie manuell 
fliegen, desto mehr scheinen sie kostspielige Zusatzschulungen zu benötigen, um verlorene 
Fähigkeiten aufzufrischen. Das Risiko solcher Ereignisse ist mit etwa einem tödlichen Unfall 
pro 2 Millionen Flügen sehr gering. 

Eine Überlegung ist, Piloten die Möglichkeit zu geben, manuelle Flugpraxis ohne Kosten für 
die Fluggesellschaften zu erhalten, indem ihnen erlaubt wird, wie IFALPA empfiehlt, die 
Kontrolle während kommerziellen Flügen zu übernehmen. Europa lehnt dies ab, denn wenn 
Fluggesellschaften es zulassen, wissen sie aus der Überwachung der Betriebsflugdaten, dass 
es mehr energiereiche oder instabile und verpasste Anflüge gibt. Die USA hingegen befür¬ 
worten manuelle Flugpraxis während kommerziellen Flügen, denn es ist sinnvoll, Piloten 
Praxis bei der manuellen Handhabung zu geben, wenn die Bedingungen sicher sind. Aller¬ 
dings ist das Wetter in den USA hinsichtlich der Flugbedingungen in der Regel gutmütiger 
und weniger schwankend als in Europa, da ein großer Teil der USA in mediterranen Breiten¬ 
graden liegt, während nur Südeuropa die meiste Zeit des Jahres gutes Wetter genießt. Dies 
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führt zu einem US-Flugverkehrsmanagementsystem, das eher bereit ist, Flüge nach Sichtflug¬ 
regeln zuzulassen. Zudem hat der größte Teil des Luftraums der USA, mit Ausnahme einiger 
wichtiger Terminalbereiche, weniger intensiven Verkehr als im dicht besiedelten Europa. 

Die FAA-Statistiken zeigen, dass in der allgemeinen Luftfahrt, dem Umfeld, in dem alle zivi¬ 
len Flieger, einschließlich der Piloten von Fluggesellschaften, ihr berufliches Lernen begin¬ 
nen, für den Kontrollverlust, in der Regel ein Strömungsabriss als die häufigste Ursache für 
tödliche Unfälle ausgemacht werden kann. 

Die Besatzungen müssen darauf vorbereitet sein, wenn die Bedingungen im Cockpit deutlich 
von den gewohnten Bedingungen abweichen. Da die im Anfangsstadium der Pilote nk arriere 
gewonnenen Erkenntnisse eher langfristig erhalten bleiben, als jedes andere erlernte Wissen, 
ist es wichtig, dass Piloten nicht nur Tests bestehen, sondern auch vollständig verstehen, was 
sie lernen, damit Ihnen Fehleinschätzungen sich nicht durch eine gesamte Karriere folgen. 
(Learmount 2017) 


Wenngleich die absolute Anzahl der fatalen Flugzeugunglücke zwar im Verlauf der letzten 
Jahrzehnte abgenommen hat (Bild 1.1), so ist prozentual häufig ein Pilotenfehler ausschlag¬ 
gebend (Tabelle 1.1). 


Tabelle 1.1 Tödliche Flugzeugunfälle nach Jahrzehnten - nach Kebabjian 2019 


Dccade 

1950$ 

1960$ 

1970$ 

19BO$ 

1990$ 

2000$ 

2010$ 

All 

Pilot Error 

50% 

53% 

49% 

42% 

49% 

50% 

57% 

49% 

Mechanical 

26% 

27% 

19% 

22% 

22% 

23% 

21% 

23% 

Weather 

15% 

7% 

10% 

14% 

7% 

8% 

10% 

10% 

Sabotage 

4% 

4% 

9% 

12% 

8% 

9% 

a% 

8% 

Othor 

5% 

9% 

13% 

10% 

14% 

10% 

4% 

10% 
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Im Flugzeug werden die Geschwindigkeiten über den mit sogenannten Pitot-Static-Systemen 
gemessenen statischen sowie totalen Druck ermittelt (Bild 1.1). Die Differenz zwischen dem 
statischen Druck und absoluten Gesamtdruck entspricht dem dynamischen Druck. Die folgen¬ 
den Zusammenhänge gelten für Unterschallgeschwindigkeiten bei M < 0,4 nach Scheiderer 
2008. 


p = p s + q 
P 2 

<7 = 2 * 

, P 2 
P = Ps + 2 V 


v 


■ (P - Ps) 


( 1 . 1 ) 

( 1 . 2 ) 

(1.3) 

(1.4) 


p Gesamtdruck 

p s statischer Druck 

q dynamischer Druck, inkompressibel 

p Dichte 

v Geschwindigkeit 
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ALTERNATE STATIC PORT 



Bild 1.2 Pitot-Static-System des CRJ 900 (Scheiderer 2008) 

Aus dem dynamischen Druck ergibt sich die angezeigte Geschwindigkeit (IAS), also die Ge¬ 
schwindigkeit, die dem Piloten auf dem Fahrmesser angezeigt wird (Bild 1.2). Die angezeigte 
Geschwindigkeit ist jedoch aufgrund von Statikdruckquellenfehlern, Ausrichtungsfehlem, 
Dichteänderungen mit der Höhe und Energieunterschieden am Flugzeugrumpf wegen Strö¬ 
mungsvorgängen fehlerbehaftet (Bild 1.3). Daher wird der tatsächliche Druck nicht exakt er¬ 
fasst. Werden die Positionierungsfehler berücksichtigt, so erhalt man die kalibrierte Ge¬ 
schwindigkeit (CAS). 


Membran 



Statikdruckleitung 


Bild 1.3 Aufbau eines Fahrtmesser (Sachs 2014) 

Werden zudem Kompressibilitätseffekte berücksichtigt, dann resultiert die äquivalente Ge¬ 
schwindigkeit (EAS). 

Wenn darüber hinaus die sinkende Dichte mit steigender Höhe Beachtung findet, so erhält 
man die Wahre Fluggeschwindigkeit (TAS). 

Neben den genannten Geschwindigkeiten, ist oft auch die Machzahl (bei hohen Geschwindig¬ 
keiten) entscheidend. Welche Fluggeschwindigkeit für den Piloten gerade relevant ist, hängt 
von dem jeweiligen Flugabschnitt, der Situation und der Fragestellung ab, beispielsweise ob 
bei einem Steigflug der Abschnitt mit konstanter angezeigter Geschwindigkeit oder mit kon- 
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stanter Machzahl geflogen wird. Im Flughandbuch sind zudem i.d.R. die Grenzwerte für die 
Geschwindigkeit in IAS angegeben, weshalb grundsätzlich mit IAS geflogen wird, da die 
Grenzwerte dann unabhängig von der Höhe gelten, sofern man sich an der angezeigten Ge¬ 
schwindigkeit orientiert. 

Für die Reichweite oder Flugzeit ist wiederum die kinematische, bzw. die Bahngeschwindig¬ 
keit (v K ) entscheidend, welche aus aerodynamischer Geschwindigkeit (v a ) und der Windge¬ 
schwindigkeit (v w ) resultiert. Nach Nguewo 2018: 

v k = v w + v a (1.5) 

v a aerodynamische Geschwindigkeit, entspricht TAS 

v k kinematische Geschwindigkeit, (Bahngeschwindigkeit) 

v w Windgeschwindigkeit 

Eine mögliche Faustformel, um die kalibrierte Geschwindigkeit in die wahre Geschwindigkeit 
umzurechnen, ergibt sich indem die Flugfläche durch 10 dividiert und das Resultat mit 6 mul¬ 
tipliziert wird. Anschließend werden der CAS-Wert sowie der TAT-Wert hinzuaddiert. 

Hierzu ein Beispiel: 

FL: 350 
GAS: 280 kt 
TAT: -15 °C 

v R0T = ■ 6 + 280 - 15) kt = 475 kt 

CAS kalibrierte Geschwindigkeit (Calibrated Airspeed) 

FL Flugfläche (Flight Level) 

TAT absolute Temperatur (Total Airtemperature) 

v ROr wahre Geschwindigkeit aus Faustformel (True Airspeed) 

Die Faustformel entstammt dem Skript von Prof. Dr. Scholz (Scholz 2019), dass an der HAW 
Hamburg ergänzend zur Vorlesung Flugmechanik bereitgestellt wird. 

Somit kann die Faustfonnel wie in Gl. (1.6) dargestellt werden. Die Faustformel aus Gl. (1.6) 
gilt es im weiteren Verlauf zu überprüfen: 
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Vrot — ~~ + v c + T t in kt (1.6) 

v c Kalibrierte Fluggeschwindigkeit (Calibrated Air Speed, CAS) in kt 

F l Flugfläche (Flight Level) 

T t absolute Temperatur (Total Air Temperature) in °C 

v rot Wahre Fluggeschwindigkeit (True Air Speed, TAS) in kt 
berechnet aus der Rule of Thumb (ROT) 


Der Index ROT (rule of thumb) wird für alle Ergebnisse aus der Faustformel verwendet. 


Korrektur von Statik- 
druckquellen-messfehler 


Korrektur flir 
Dichte 



Bild 1.4 Geschwindigkeitsabhängigkeiten (nach Scheiderer 2008) 


1.2 Begriffsdefinitionen 

Eine Faustformel ist „die grobe einfache Formel, mit der man eine überschlägige Berech¬ 
nung durchführen kann.“ (Scholze-Stubenrecht 2015) 

Als synonymer Begriff für Faustformel kann aus dem englischen abgeleitet Daumenregel 
(rule of thumbs) verwendet werden. 

Piloten benötigen oft keine exakten Angaben, sondern vielmehr Richtwerte mit gewisser To¬ 
leranz, die sich auf die abgelesenen Werte beziehen. Dazu werden dann oft entsprechende 
Faustformeln verwendet, die ausreichend genaue Ergebnisse liefern und dabei vor allem den 
Vorteil ergeben schnell eine Lösung anzubieten und gleichzeitig einen überschaubaren Kopf¬ 
rechenaufwand fordern. 
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Faustformeln gibt es in der Luftfahrt für eine Vielzahl von Parametern, so z.B. um u.a. ver¬ 
schiede Höhen, Einheiten oder Geschwindigkeiten umzurechnen. 

Die Flugmechanik befasst sich mit der Anwendung der Gesetze der Mechanik und Dynamik 
auf Flugzeuge. Aufgabe der Flugmechanik ist es, die am Flugzeug wirkenden Kräfte und 
Momente zu bestimmen, und daraus den Flugzustand mit seinen Zustandsgrößen und die Sta¬ 
bilität des Flugzeugs zu berechnen. (Klußmann 2007) 


1.3 Ziel der Arbeit 


Ziel dieser Projektarbeit ist es, die Kopfrechenmethode zur Umrechnung von CAS in TAS aus 
Gl. (1.6) auf ihre Gültigkeit zu prüfen, respektive einen Gültigkeitsbereich zu definieren. 


v rot — 


6 -_Fl 
10 


+ v c + T t 


(1.6) 


Hierzu ist es notwendig eine hinreichend hohe An z ahl an TAS-Werten über die Faustfonnel 
aus den zu prüfenden kalibrierten Geschwindigkeiten und Höhen zu ermitteln und mit den be¬ 
reits bewährten Gleichungen der Flugmechanik mittels Zusammenhangs zwischen Tempera¬ 
tur, Höhe, CAS, TAT, Machzahl und TAS bestenfalls zu verifizieren oder ggfs, zu widerle¬ 
gen. 


1.4 Literaturübersicht 

Die Definitionen beziehen sich zumeist auf das Lexikon der Luftfahrt“ von Klußmann 2007. 
Dort stehen die meisten für diese Studienarbeit relevanten Fachtermini der Luftfahrt. 

Die wichtigste und die an der häufigsten genutzten Quelle, insbesondere in Bezug auf die 
verwendeten Gleichungen, ist das Manuskript „Flight Mechanics Lecture not es“ von Y- 
oung 2001. Das Skript wird begleitend für die Vorlesung „Flugmechanik 1“ an der HAW von 
Prof. Scholz genutzt und bietet einen sehr guten Überblick über die wichtigsten flugmechani¬ 
schen Zusammenhänge und Gleichungen. 

Ergänzend zum Manuskript „Flight Mechanics Lecture notes “ von Young 2001 werden „An¬ 
gewandte Flugleistungen “ von Scheiderer 2008, „ Handbuch der Luftfahrzeugtechnik “ von 
Sachs 2014 sowie das Vorlesungsskript. „ Vorlesung Flugmechanik 1 “ von Nguewo 2018 
verwendet. Weitere Quellen finden teils nur eimnalig Verwendung und werden an dieser Stel¬ 
le nicht weiter aufgeführt, finden sich jedoch im Literaturverzeichnis. 
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1.5 Aufbau der Arbeit 

Um diese Projektarbeit in vollem Umfang nachvollziehen zu können, ist es hilfreich, wenn¬ 
gleich nicht zwingend erforderlich, die interaktiven Diagramme und Tabellen, sowie das Li- 

ve-Skript aus Matlab zu nutzen. Hierbei werden Vorkenntnisse der Programme jeweils vo¬ 
rausgesetzt. 

Im Hauptteil der Projektarbeit steht die Verifizierung der Faustformel aus Gl. (1.6) im Fokus. 

Abschnitt 2 beschreibt die grundsätzlichen aus flugmechanischem Aspekt relevanten 
Zusammenhänge und Parameter der Erdatmosphäre und die daraus resultie¬ 
renden Gleichungen. 

Abschnitt 3 erklärt den Aufbau des „Comressibility Correction Chart“ und es werden an 
Anwendungsbeispielen exemplarisch einige CAS-Werte in TAS-Werte 
überführt. 

Abschnitt 4 stellt die zu Grunde liegenden Gleichungen und die Beziehungen zu den ge¬ 
nerierten Diagrammen in Excel dar. Es wird ein Zusammenhang zwischen 
den bewährten Gleichungen der Flugmechanik und der Faustfonnel aus Gl 
(1.6) sowie dem daraus ermittelten relativen Fehler der Faustformel in Pro¬ 
zent hergeleitet. 

Abschnitt 5 zeigt die Zusammenhänge der nachgebildeten ISA-Tabellen und den Tabel¬ 
len für den relativen Fehler der Faustformel in Abhängigkeit von Höhe und 
Geschwindigkeit. 

Abschnitt 6 befasst sich mit der Programmierung eines Live-Skripts in Matlab zur Visu¬ 
alisierung von zwei -und dreidimensionalen Diagrammen in Abhängigkeit 
von den Parametern CAS, H und A T um die relativen Abweichungen der 
Faustformel darstellen zu können. 

Abschnitt 7 präsentiert die Ergebnisse aus den Tabellen, der Diagramme und aus dem 
Matlab-Skript. 



20 


2 Internationale Standard Atmosphäre (ISA) 


2.1 Erdatmosphäre 


Das Leistungsvermögen eines Flugzeuges ist im hohen Maß abhängig von den vorherrschen¬ 
den atmosphärischen Bedingungen. Die Atmosphäre wird durch Sonneneinstrahlung erhitzt. 
Die Temperatur der bodennahen Luftschichten erhöht sich dadurch. Die erwärmte Luft steigt 
auf und kühlt sich dabei ab. Ist keine weitere Abkühlung mehr möglich, so erreicht man die 
Tropopause. Die Erdatmosphäre lässt sich über den Temperaturverlauf in vier Schichten auf¬ 
teilen (Bild 2.1). Die unterste Schicht der Erdatmosphäre bezeichnet die Troposphäre. Dar¬ 
über liegt die Stratosphäre, wobei die Tropopause die Grenzschicht zwischen der Troposphäre 
und der Stratosphäre bildet. In der Stratosphäre ist die Temperatur weitestgehend konstant. 
Die Tropopausenhöhe variiert in den unterschiedlichen Breiten und ist aufgrund ihrer Grenz- 
schichtfünktion für die Flugplanung und Flugdurchführung von besonderer Bedeutung. Die 
reale Atmosphäre unterliegt ständigen Änderungen. Wechselnde Jahreszeiten, tatsächlich 
vorherrschende Winde und Luftmassenbewegungen haben großen Einfluss auf die Tempera¬ 
tur, den Luftdruck und die Dichte. (Scheiderer 2008) 




-80 -60 -40 -20 0 +20 T/°C 


Bild 2.1 Temperaturverlauf der Atmosphäre (Demtröder 1994, S. 207) 


Zwecks Vergleichbarkeit der Flugleistungs- und Triebwerksparameter ist eine einheitliche 
Basis notwendig. Zudem erfordert die (Höhen-) Staffelung der Flugzeuge durch die Flugver¬ 
kehrskontrolle einheitliche Referenzwerte bzgl. Druck, Dichte und Temperatur, damit diese 
die gleichen Werte angezeigt bekommen. Die Internationale Standard Atmosphäre (ISA) der 
ICAO bietet eine solche Referenz. 
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2.2 Allgemeine Zusammenhänge der Atmosphäre 


Zustandsgleichung: 


V 

p = Rl ' t 


Po = 


Po 


R l -t 0 


Hydrostatische Gleichung: 


dp 

dh 


-p-g 


Aus Gl. (2.1) und Gl. (2.3) ergibt sich 


rP 1 r 

- dp = - 
KoP Jh 


g 


)p0 P J HO Rl ‘ T 


dH 


( 2 . 1 ) 

( 2 . 2 ) 


(2.3) 


(2.4) 


g Erdbeschleunigung 

H Geopotenzielle Höhe 

p 0 Referenzdruck auf Meereshöhe 

R l Gaskonstante, Luft 

T Temperatur 

T {) Referenztemperatur auf Meereshöhe 

Während in der Troposphäre ein konstanter Temperaturfall von 6,5 °C pro km bei steigender 
Höhe vorausgesetzt wird, besteht in der Stratosphäre annähernd Isothermie bis zu einer Höhe 
von 20 km. Daher muss die Integration von Gl. (2.4) jeweils für beide Schichten der Atmo¬ 
sphäre getrennt ausgeführt werden. 


2.3 Zusammenhänge für die Troposphäre 


Die Integration von Gl. (2.4) ergibt bei konstantem Temperaturgradienten L: 


8 


P_ 

Po 





g 

RL 


(2.5) 
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Unter Beachtung der tatsächlichen vorherrschenden Temperaturverhältnisse ergibt sich eine 
Differenz zu der Referenztemperatur, die durch AT ausgedrückt werden kann, woraus Gl. 
(2.6) resultiert: 


o_ 


P _( L-H y-L 

(2.6) 

Po V T 0 + AT) 

Außerdem gilt: 


5 


o - — 

(2.7) 

0 

8 Druckverhältnis 

9 Temperaturverhältnis 

er Dichteverhältnis 


AT Differenz aus Referenztemperatur und tatsächlicher Temperatur auf Meereshöhe 


2.4 

Zusammenhänge für die Stratosphäre 


Die Integration von Gl. (2.4) ergibt bei konstanter Temperatur T: 



— — e r ' t t < ' H Ht '•* 

Pt 

(2.8) 

H r 

Höhe der Tropopause 


Pt 

Referenzdruck in Höhe der Tropopause 


T t 

Referenztemperatur in Höhe der Tropopause 


Für die Tropopause gelten zudem folgend aufgeführte Verhältnisse: 



T t 

0 — — — 0,75187 

T To 

(2.9) 


5 =?L = 0,22336 

Po 

(2.10) 


a T = — = 0,29708 

Po 

(2.11) 
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Außerdem gilt: 



( 2 . 12 ) 


8 t Druckverhältnis, Tropopause 

Q t Temperaturverhältnis, Tropopause 

o T Dichteverhältnis, Tropopause 

p T Referenzdichte in Höhe der Tropopause 

Po Referenzdichte auf Meereshöhe (MSL) 

Anmerkung: Die angegebenen Werte aus Gl. (2.9) bis Gl. (2.11) gelten nur für den Fall, dass 
ideale ISA-Bedingungen vorliegen, folglich A7=0 ist. Für abweichende Temperaturverhält¬ 
nisse werden auch in der Stratosphäre im weiteren Verlauf die Grenzwerte entsprechend an¬ 
gepasst. 


7> = -56,5 + A T 


(2.13) 


2.5 Konstanten in der Atmosphäre 

Die im Konstanten müssen zunächst in Non-SI Einheiten umgerechnet werden. Längenanga¬ 
ben erfolgen hierbei in ft, Geschwindigkeitsangaben in kt. 


lm = 3,280839895 ft 

(2.14) 

1ms“ 1 = 1,94384494 kt 

(2.15) 

ms~ 2 = 3,280839895 fts~ 2 

(2.16) 


Tabelle 2.1 Allgemeine Konstanten 


Bezeichnung 

Symbol 

Wert 

SI-Einheit 

Wert 

Non-SI 

d. physikalischen Größe 





Einheit 

Isentropenexponent (Luft) 

K 

1,4 


1,4 


Spez. Gaskonstante (Luft) 

Rl 

287,053 

K -1 m 2 s“ 2 

3089,813 

K _1 ft 2 s“ 2 

Gravitationskonstante 

go 

9,80665 

ms“ 2 

32,174 

CN 

1 

d=5 

Erdradius 

Tearth 

6371■10 3 

m 

20,9 ■ 10 6 

ft 
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Tabelle 2.2 Konstante Parameter der Troposphäre nach ISA 


Bezeichnung 

Symbol 

Wert 

SI- 

Wert 

Non-SI 

d. physikalischen Größe 



Einheit 


Einheit 

(Referenz-) Temperatur (MSL) 

T 0 

288,15 

K 

- 


(Referenz-) Temperatur (MSL) 

To 

15 

C° 

- 


T emperaturgradient 

L 

6,5 ■ IO" 3 

Km -1 

1,9812 ■ 10 -3 

K ft -1 

Schallgeschwindigkeit (MSL) 

a 0 

340,294 

mD' 1 

661,48 

kt 


Tabelle 2.3 Konstante Parameter der Stratosphäre nach ISA 


Bezeichnung 

Symbol 

Wert 

SI- 

Wert 

Non-SI 

d. physikalischen Größe 



Einheit 


Einheit 

(Referenz-) Temperatur (Tropopause) 

Ts 

216,65 

K 

- 


(Referenz-) Temperatur (Tropopause) 

Ts 

-56,5 

C° 

- 


Temperaturgradient (bis ca. 20 km) 

Ls 

0 

Km -1 

0 

Kft -1 

Tropopausenhöhe 

H T 

11000 

m 

36089,24 

ft 
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3 Comressibility Correction Chart 

3.1 Einführung in das Chart 


Ziel dieser Arbeit ist es, die Faustformel zur Umrechnung von CAS in TAS zu verifizieren, 
bzw. einen Gültigkeitsbereich zu definieren. Hierzu muss jedoch zunächst eine Methode ge¬ 
wählt werden, um die zu prüfenden Werte der Faustformel zu bestätigen, respektive um die 
Abweichung quantitativ zu benennen. Eine solche Vergleichsmöglichkeit bietet das „com- 
pressibility correction chart“ aus Bild 3.1. Bei Geschwindigkeiten ab ca. Mach 0,3 darf die 
Luft nichtmehr wie für Gl. (3.2) als inkompressibel angenommen werden. Für inkompressible 
Fluide ist stattdessen Gl. (3.3) für den dynamische Druck zu verwenden. Im Chart wird für ei¬ 
nen CAS-Wert die Differenz zu dem jeweiligen EAS-Wert in Abhängigkeit von der Höhe als 
AV in kt abgelesen. Bei den Höhe nk urven handelt es sich um die Drukchöhe h P , wobei unter 
ISA-Bedingung nach Gl. (3.1) für ZF = 0 gilt, dass h p — H. 


hp _ T 0 

77 ~~ T 0 + AT 


(3.1) 



COMPRESSIBILITY CORRECTION 


-30 


O 
lii 
GC 

Sog -20 

I I I v V 


LU 

Q. 

cn 

GC 

< 


2 

* 

ö 

> 

< 


-10 


0 

100 


200 300 400 500 

CALIBRATED AIRSPEED, KNOTS 


600 


Bild 3.1 


Compressibility Correction Chart (Dole 2017) 
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3.2 Compressibility Correction Chart in Excel 

Um exemplarisch einige Beispiele zu visualisieren wird zunächst das Diagramm in Excel 
übertragen. Hierzu werden die Zusammenhänge aus Gl. (3.2) bis Gl. (3.9) genutzt, um die 
Höhe nk urven mit einem Intervall von 5000 ft zwischen 5000 ft und 65000 ft und Machzahl¬ 
kurven im Intervall von 0,1 Mach im Bereich 0,4 Mach bis 1 Mach darzustellen. Das Chart ist 
interaktiv und zeigt inkl. Hilfslinien über den gewählten CAS-Wert sowie der gewählten Hö¬ 
he den gesuchten AU-Wert an. Das Originalchart stellt negative AU-Werte in Abhängigkeit 
von CAS dar. In dem interaktiven Chart in Excel hingegen werden zwecks einfacherer Hand¬ 
habung die Werte positiv abgelesen, wodurch sich dadurch Gl. (3.8) ergibt. Würde hingegen 
mit dem Originalchart aus Bild 3,1 gearbeitet, so würde der abgelesene AU-Wert addiert, statt 
subtrahiert. 


R 



(3.2) 


Rc = V 


1) p(//) 

2 ■ k p(H) 


i K-l 


TAS 2 + 1 


- 1 


(3.3) 


TAS = 



(3.4) 


AU =| EAS- CAS | 


(3.5) 


AV=!M- a - Vö-CAS \ 


(3.6) 
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M = 


5 < 


4 


1 + 


y — 1 (v , 


a 0 / 


lA 1 

1 

Y~ 1 

Y 

- 1 

• +1 

- 1 

J 

■ 

j 


(3.7) 


£715 = CAS - AE 


(3.8) 


EAS 
TAS = —= 
yjo 


(3.9) 


3.3 Anwendungsbeispiele aus dem Chart 

Es werden folgend exemplarisch 3 Beispiele unter Anwendung des “Compressibility Correc- 
tion Chart“ aufgezeigt. Als Input Parameter dienen der CAS-Wert sowie die Druckhöhe, die 
unter der Annahme, dass ISA-Bedingungen herrschen der geopotenziellen Höhe entspricht. 
Eine kurze Beschreibung zum Aufbau des Charts befindet sich im Anhang auf die bei Nut¬ 
zung der Originaldatei zurückgegriffen werden kann. Zu sehen sind in Blau die Höhe nk urven 
und in Rot die Machzahlkurven. Auf der vertikalen Achse können jeweils die AE-Werte abge¬ 
lesen werden. 
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compressibility correction to CAS 



Bild 3.2 Compressibility Correction Chart, Beispiel 1 

Es wird mit einem CAS-Wert von 280 kt und einer Höhe von 40000 ft aus Bild 3.2 ein AV- 
Wert von 15,82 angezeigt. Dieser wird auf 16 gerundet. Zunächst wird mittels Gl. (3.8) EAS 
berechnet. Darauffolgend wird der entsprechende o -Wert aus der ISA-Tabelle für 40000 ft 
mit 0,24617 abgelesen (Tabelle 3.1) und schlussendlich aus Gl. (3.9) der TAS-Wert ermittelt. 
Bei den Beispielen 2 und 3 aus Bild 3.3 sowie Bild 3.4 wird analog verfahren. 


Tabelle 3.1 Auszug ISA-Tabelle, 40000 ft 


Höhe 

FL 

P 

P 

T 

a 

6 

G 

CD 

[ft] 

[km] 

[ft/100] 

[hPa] 

[kg/m 3 ] 

[°C] 

[K] 

[kt] 




40.000 

12,19 

400,00 

187,539 

0,3016 

- 56,5 

216,6 

573,57 

0,18509 

0,24617 

0,75187 


Beispiel 1: 

CAS = 250 kt 
H = 40000 ft 

er 40 ooo = 0,24617 (Tabelle 3.1) 
ÄV ~ 16 kt 

EAS = (250 - 16) kt = 234 kt 
234 

TAS = . kt » 472 kt 

V0,24617 
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AV 


compressibility correction to CAS 


35,00 


30,00 


25,00 


20,00 


15,00 


10,00 


5,00 


0,00 



CAS [kt] 


Bild 3.3 Compressibility Correction Chart, Beispiel 2 


Tabelle 3.2 Auszug ISA-Tabelle, 35000 ft 


Höhe 

FL 

P 

Q. 

T 

a 

5 

G 

e 

[ft] 

[km] 

[ft/100] 

[hPa] 

[kg/m 3 ] 

[°C] 

[K] 

[kt] 




35.000 

10,67 

350,00 

238,42 

0,3796 

- 54,3 

218,8 

576,42 

0,23531 

0,30988 

0,75935 


Beispiel 2: 

CAS = 300 kt 
H = 35000 ft 

<735000 — 0,30988 (Tabelle 3.2) 
AV * 20 kt 

EAS = (300 - 20) kt = 280 kt 
280 

TAS ~ — kt ~ 503 kt 

V0,30988 
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AV 

40,00 

35,00 


compressibility correction to CAS 


25,00 


15,00 



Bild 3.4 


CAS [kt] 

Compressibility Correction Chart, Beispiel 3 


Beispiel 3: 

CAS = 280 kt 
H = 35000 ft 

0-35000 — 0,30988 (Tabelle 3.2) 

AV * 16,5 kt 

EAS = (280 - 16,5) kt = 263,5 kt 
263,5 

TAS = — kt « 473 kt 

V0,30988 


Um die Faustformel mithilfe des Charts zu verifizieren, müssten jeweils eine Vielzahl an 
Werten nach vorangegangenem Schema mittels Charts mit den Ergebnissen aus der Faustfor¬ 
mel abgeglichen werden. Um einen gesamten Gültigkeitsbereich zu definieren und die Ergeb¬ 
nisse später besser visualisieren zu können, wird in den folgenden Kapiteln eine weitere Fler- 
angehensweise erarbeitet. 
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4 Herleitung des relativen Fehlers 

Zunächst wird ein direkter mathematischer Zusammenhang zwischen der Faustformel und 
den bewährten Gleichungen der Flugmechanik hergestellt. Hierbei ist es notwendig, die 
Faustformel und die flugmechanisch relevanten Gleichungen in Abhängigkeit derselben vari¬ 
ablen Parameter auszudrücken. 

Es wird davon ausgegangen, dass die Geschwindigkeit in Form der kalibrierten Geschwindig¬ 
keit bereits bekannt ist. Nachfolgend aufgeführte Gleichungen, sofern nicht explizit anderwei¬ 
tig genannt, entstammen Young 2001. 

Die Machzahl berechnet sich wie in Gl. (3.7) angegeben: 


ao 

M 

vc 

y 


M = 


5 < 




1 + 


y — 1 (v t 


Vc\ 

aj 


lA 1 

1 

Y-l \ 

Y 

- 1 

• +1 

- 1 

J 

■ 

j 


Schallgeschwindigkeit auf Meereshöhe 
Machzahl 

kalibrierte Geschwindigkeit (calibrated airspeed, CAS) 
Wärmekapazitätsverhältnis 


(3.7) 


Das Druckverhältnis lässt sich in Abhängigkeit der Temperatur folgend darstellen: 


ö = 



L-H 
T 0 + ATJ 


(4.1) 


Anmerkung: Es wird im weiteren Verlauf ausschließlich mit der geopotenziellen Höhe H ge¬ 
rechnet. Diese hat den Vorteil, dass g als konstant angenommen werden darf. Der Zusam¬ 
menhang der geometrischen Höhe h und der geopotenziellen Höhe H wird in Gl (4.2) aufge¬ 
führt: 


h = 


' earth 


H 


f earth H 


(4.2) 


Tearth 

h 


Erdradius 

Geometrische Höhe 
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Die wahre Geschwindigkeit ergibt sich aus der Machzahl und der Schallgeschwindigkeit: 

v — a - M (4.3) 

v wahre Geschwindigkeit (true airspeed, TAS) 

a Schallgeschwindigkeit 

Sowohl die Schallgeschwindigkeit als auch die Machzahl sind Fu nk tionen in Abhängigkeit 
von der Höhe, respektive von der Temperatur und des Drucks in der jeweiligen (Flug-) Höhe. 

Allgemein gilt für die Schallgeschwindigkeit: 

a(H) = yjT(H) ■ R l ■ Y (4.4) 


T(H) Temperatur in Höhe H 

a(H) Schallgeschwindigkeit in Höhe H 


Stellt man das Verhältnis zu der Schallgeschwindigkeit auf Meereshöhe dar, so ergibt sich: 


a(H) _ y/T(H~) ■ R l -Y 
a o JT 0 -R l -Y 

a 0 Schallgeschwindigkeit auf Meereshöhe 

Daraus folgt durch Kürzen: 


a(H ) - a 0 



Die Temperatur ergibt sich aus: 


T(H) = T 0 + AT — L ■ H 

Werden Gl. (4.6) und Gl. (4.7) zusammengefasst, so resultiert: 


(4.5) 


(4.6) 


(4.7) 





a(H ) = a 0 


(4.8) 
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Werden nunmehr Gl. (3.7) und Gl. (4.8) in Gl. (4.3) überführt, so erhält man die gesuchte 
wahre Geschwindigkeit in Abhängigkeit von der Höhe H, der Temperaturdifferenz A T und 
der kalibrierten Geschwindigkeit v c : 


, fAT — L ■ H\ 

v = a ° l 1 + ( r„ ) 


5 < 


4 


1 + 


y — 1 (V 


W\ 

aj 




Y~ 1 \ 

DA 3 

■ 

Y 

y - 1 



~ i 

• + 1 

- 1 

J 

■ 



(4.9) 


Das Druckverhältnis ö berechnet sich nach Gl. (4.1). 


Nun gilt es die Faustfonnel über die gleichen variablen Parameter darzustellen. 

6 

v = Fl ' w + Vc + Tt 

Damit die Ergebnisse später direkt in den Tabellen und Diagrammen vergleichbar werden und 
entsprechend über die gleichen Parameter, bzw. (Diagramm-) Achsen dargestellt werden kön¬ 
nen, müssen die Flugfläche und die wahre Temperatur ebenfalls über die Höhe H, der Tempe¬ 
raturdifferenz und der kalibrierten Geschwindigkeit ausgedrückt werden. 


Die Flugfläche ergibt sich abgeleitet aus Gl. (3.1) und Gl. (4.11) wie folgt: 


FL 


hp 

1ÖÖ 


(4.10) 


FL 


H ' 7p 

100 (T 0 + A T) 


(4.11) 


Die absolute Temperatur berechnet sich aus: 


TAT = (T 0 - L ■ H + AT ) (1 + k r 



(4.12 

2 ) 


k r Recovery Faktor 
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Werden Gl. (1.6), Gl. (4.11) und Gl. (4.12) zu einer Gleichung zusammengefasst, so resultiert 
die gesuchte Gleichung Gl. (4.13) in Abhängigkeit der variablen Parameter H, A T und v c . 


TASt/nr — 


H ‘ T 0 


R0T 100 ■ (T 0 + Ar) 10 

l f r 

Y — 1 1 

1 + V 5 5 


— + v c + (T„-L-H + AT) 




1 + 


Y — 1 (v 


aj 


y \ i 

y-1 

\ 


y 


-1+1 

- 1 

>■ 

J 

y 

/ 


(4.13) 


Das Druckverhältnis 5 berechnet sich nach Gl. (4.1). Es wird für Gl. (4.12) vorausgesetzt, 
dass kr= 1. 


Verwendete Konstanten für Luft in der Troposphäre: 

Y = 1,4 

g 0 = 9,80665 m s -1 

L = 6,5 K/km = 0,0019812 K/ft 

T 0 = 288,15 K 

R l = 287,053 J/(kgK) 

a 0 — 661,48 kt 

Des Weiteren können einige Therme zusammengefasst werden: 


y-l 

2 =°’ 2 

(4.14) 

—t— = 3,5 
y-l 

(4.15) 

y — 1 

-5 = 1 

2 

(4.16) 

y — 1 1 

y “ 3,5 

(4.17) 

9 = 5,25588 

L ' L) 

(4.18) 
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Werden Gl. (4.9) sowie Gl. (4.13) nun unter Beachtung der Konstanten Werte und Gl. (4.14) 
bis Gl. (4.18) dargestellt, so resultieren daraus 


v — 661,48 1 + 


A7 - 0,0019812 /G 


288,15 



(4.19) 


und 


v ROT 


288,15 H 


100 (288,15 + Ar) 10 


— + v c + (288,15 - 1,9812 ■ 10~ 3 H + AT ) 



(4.19) 


mit jeweils dem Druckverhältnis aus Gl. (4.1) 


ö - 



0,0019812 tf^ 5 ' 2558 
288,15 + ATJ 


(4.20) 


Nun stehen durch Gl. (4.19) und Gl. (4.20) zwei Gleichungen in Abhängigkeit derselben Pa¬ 
rameter zur Verfügung. Diese ennöglichen es über die Werte der ermittelten wahren Ge¬ 
schwindigkeit aus der Faustformel und den errechneten wahren Geschwindigkeiten aus den 
Gleichungen der Flugmechanik einen relativen Fehler in Prozent zu bestimmen. 


v - TAS , 


ROT 


v 


100 


(4.21) 


e relativer Fehler in [%] 

Zu beachten ist sowohl bei der Zusammenfassung der Konstanten als auch beim Verwenden 
der daraus resultierenden Gleichungen, dass nicht durchgängig SI-Einheiten verwendet wer¬ 
den und alle Werte einheitenkonform einzusetzen sind. Entfernungen, bzw. Höhen erfolgen in 
ft, Geschwindigkeitsangaben in kt und Temperaturen werden in Kelvin (K) angegeben. 
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5 Relativer Fehler in Excel 

5.1 Integrieren der ISA-Tabellen in Excel 

Wenngleich es fertige Auszüge mit den für diese Projektarbeit relevanten Parameter der ISA- 
Tabellen gibt, wie z.B. in Tabelle 5.1 exemplarisch dargestellt, so wäre es äußerst mühsam, 
diese manuell zu übertragen. Über die vorangegangenen Gleichungen und definierten Kon¬ 
stanten aus Kapitel 2 bis 4 kann durch übertragen der Gleichungen unter Beachtung der defi¬ 
nierten Konstanten auch ein Auszug in Excel generiert werden, der alle für die weiteren Be¬ 
rechnung relevanten Parameter der ISA-Tabelle enthält. Das Ergebnis wird in Tabelle 5.2 
gezeigt. 
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Tabelle 5.1 Auszug einer ISA-Tabelle (Scheiderer 2008, S. 101) 


Höhe (h) 

Temperatur (T) 

Druck |p| 

Dichte (p) 

ft 

m 

K 

•F 

•c 

0 


in Hg 

hPa 

6 

kg/m 3 

o 


0 

288.16 

59.0 

15.0 

1.0000 

1.0000 

29 921 

1013.25 

1.0000 

1.2250 

1.0000 

1 000 

305 

286.18 

55.4 

13,0 

0 9931 

0.9966 

28 855 

977,17 

0.9644 

1.1896 

0.9711 

2 000 

610 

284.20 

51.9 

11.0 

0 9662 

0.9931 

27 821 

942.13 

0.9298 

1.1549 

0.9428 

3000 

914 

282.22 

48 3 

9.1 

0.9794 

0.9896 

26 816 

908.12 

0.8962 

1.1210 

0.9151 

4 000 

1.219 

280.24 

44 7 

7.1 

0.9725 

0.9862 

25 842 

875.11 

0.8637 

1.0879 

0.8681 

5000 

1 524 

278.25 

41 2 

5.1 

0 9656 

0.9827 

24 896 

843.07 

0.8320 

1.0555 

0.8617 

6000 

1 829 

276,27 

37.6 

3.1 

0.9587 

0.9792 

23978 

812.00 

0.8014 

1.0239 

0.8359 

7000 

2.134 

274,29 

34.0 

1.1 

0.9519 

0,9756 

23 068 

781.85 

0.7716 

0.9930 

0.8106 

8 000 

2 438 

272.31 

30.5 

•0.8 

0.9450 

0.9721 

22225 

752,62 

0.7428 

0.9629 

0 7860 

9 000 

2.743 

270.33 

26.9 

•2.8 

0.9381 

0.9666 

21 368 

724.28 

0.7148 

0.9334 

0.7620 

10 000 

3 048 

268.35 

23.3 

-4.8 

0.9312 

0.9650 

. 

696.82 

0.6677 

0.9046 

0,7385 

11 000 

3.353 

266.37 

19.8 

•6.8 

09244 

0.9614 

19.791 

670,20 

0.6614 

0.8766 

0,7156 

12.000 

3657 

264.39 

18.2 

•8.8 

0.9175 

0,9579 

19 029 

644.41 

0.6360 

0.8491 

0.6932 

13000 

3962 

262.40 

12.6 

-10.8 

0.9106 

0.9543 

18 292 

619.43 

0.6113 

0.8224 

0.6713 

14 000 

4267 

260.42 

9,1 

•12.7 

0.9037 

0.9507 

17 577 

595.24 

0.5875 

0.7963 

0.6500 

15 000 

4 572 

258.44 

5.5 

•14.7 

08969 

0.9470 

16 866 

571,82 

0.5643 

0,7708 

0.6292 

16000 

4 877 

256.46 

1.9 

•16.7 

0 8900 

0.9434 

16.216 

549.15 

0.5420 

0.7460 

0.6090 

17 000 

5.181 

254.48 

• 16 

•18.7 

0 8631 

0,9397 

15 569 

527,22 

0.5203 

0.7218 

0.5692 

18 000 

5486 

252.50 

•5.2 

■20.7 

0.8762 

0.9361 

14 942 

506.00 

0.4994 

0.6981 

0.5699 

19 000 

5 791 

250.52 

8 8 

•22.6 

0 8694 

0.9324 

14 336 

465.48 

0.4791 

0.6751 

0.5511 

20 000 

6 096 

248.54 

■12.3 

•24.6 

18625 

0.9267 

13.750 

465.63 

0.4595 

0.6527 

0.5328 

21 000 

6 401 

246.55 

-15.9 

•26.6 

0.8556 

0.9250 

13.184 

446.45 

0.4406 

0.6308 

0.5150 

22 000 

6 705 

244,57 

• 19.5 

-28.6 

0 8487 

0.9213 

12 636 

427,91 

0.4223 

0.6095 

04976 

23 000 

7.010 

242.59 

•23.0 

-30.6 

18419 

0,9175 

12.107 

410.01 

0.4046 

0.5888 

0.4807 

24 000 

7.315 

240.61 

■266 

■32.5 

0.8350 

0.9138 

11 597 

392.71 

0.3676 

0.5686 

0 4642 

25 000 

7620 

236.63 

-30 2 

■34.5 

0.8281 

0.9100 

11.103 

376,01 


- : - 

0 4481 

26 000 

7925 

236.65 

•33.7 

-36.5 

0 8212 

0.9062 

10627 

359.89 

0.3552 

0.5298 

0.4325 

27 000 

8229 

234.67 

•37.3 

-38.5 

0.8144 

0.9024 

10.168 

344.33 

0.3398 

0.5112 

0.4173 

28.000 

8 534 

232.69 

•40 9 

•40.5 

0 8075 

0.8966 

9.725 

329.32 

0.3250 

0.4931 

0.4025 

29 000 

8839 

230,71 

•44 4 

-42.5 

0 8006 

0.6948 

9 297 

314.85 

0.3107 

0.4754 

0,3881 

30 000 

9 144 

228,72 

-4-8 0 

-44.4 

0.7937 

0.8909 

8 865 

300.90 

0.2970 

0.4583 

.3741 

31 000 

9448 

226,74 

•51.6 

•46.4 

0 7869 

0.8870 

8 488 

287,45 

0.2837 

0,4417 

0.3605 

32 000 

9 753 

224.76 

■55.1 

-48.4 

0.7800 

0.8832 

8 106 

274.49 

0.2709 

0.4255 

0.3473 

33 000 

10058 

222.78 

-58.7 

-50.4 

0.7731 

0.8793 

7.737 

262.01 

0.2586 

0.4097 

0.3345 

34 000 

10 363 

220.80 

-622 

-52.4 

0.7662 

0,8753 

7 382 

249.99 

0.2467 

0.3944 

0.3220 

35 000 

10 668 

218.82 

-658 

•54.3 

0.7594 

0.6714 

7.041 

238.42 

- 

0.3796 

C 3099 

36 000 

10 972 

216.84 

-694 

•56.3 

0.7525 

0.8675 

6 712 

227.26 

0.2243 

C 3652 

0.2981 

36.069 

11.000 

216.66 

<9.7 

-56.5 

0.7519 

0.8671 

6.863 

226.32 

0.2234 

0.3639 

0.2971 
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Tabelle 5.2 ISA-Tabelle in Excel generiert 


H 

ft 

km 

FL 

ft/100 

P 

hPa 

P 

kg/m 3 

T 

°C 

K 

a 

kt 

6 

o 

0 

0 

0 

0 

1 . 013,25 

1,2250 

15,0 

288,15 

661,47874 

1,0000 

1,0000 

1,0000 

1.000 

0,30 

10,00 

977,17 

1,1896 

13,0 

286,2 

659,20 

0,9644 

0,9711 

0,9931 

2.000 

0,61 

20,00 

942,13 

1,1549 

11,0 

284,2 

656,91 

0,9298 

0,9428 

0,9862 

3.000 

0,91 

30,00 

908,12 

1,1210 

9,1 

282,2 

654,62 

0,8962 

0,9151 

0,9794 

4.000 

1,22 

40,00 

875,11 

1,0879 

7,1 

280,2 

652,32 

0,8637 

0,8881 

0,9725 

5.000 

1,52 

50,00 

843,07 

1,0555 

5,1 

278,2 

650,01 

0,8320 

0,8617 

0,9656 

6.000 

1,83 

60,00 

812,00 

1,0239 

3,1 

276,3 

647,69 

0,8014 

0,8359 

0,9587 

7.000 

2,13 

70,00 

781,85 

0,9930 

1,1 

274,3 

645,36 

0,7716 

0,8106 

0,9519 

8.000 

2,44 

80,00 

752,62 

0,9629 

- 0,8 

272,3 

643,03 

0,7428 

0,7860 

0,9450 

9.000 

2,74 

90,00 

724,28 

0,9334 

- 2,8 

270,3 

640,69 

0,7148 

0,7620 

0,9381 

10.000 

3,05 

100,00 

696,82 

0,9046 

- 4,8 

268,3 

638,33 

0,6877 

0,7385 

0,9312 

11.000 

3,35 

110,00 

670,20 

0,8766 

- 6,8 

266,4 

635,97 

0,6614 

0,7156 

0,9244 

12.000 

3,66 

120,00 

644,41 

0,8491 

- 8,8 

264,4 

633,60 

0,6360 

0,6932 

0,9175 

13.000 

3,96 

130,00 

619,43 

0,8224 

- 10,8 

262,4 

631,22 

0,6113 

0,6713 

0,9106 

14.000 

4,27 

140,00 

595,24 

0,7963 

- 12,7 

260,4 

628,84 

0,5875 

0,6500 

0,9037 

15.000 

4,57 

150,00 

571,82 

0,7708 

- 14,7 

258,4 

626,44 

0,5643 

0,6292 

0,8969 

16.000 

4,88 

160,00 

549,15 

0,7460 

- 16,7 

256,5 

624,03 

0,5420 

0,6090 

0,8900 

17.000 

5,18 

170,00 

527,22 

0,7218 

- 18,7 

254,5 

621,62 

0,5203 

0,5892 

0,8831 

18.000 

5,49 

180,00 

506,00 

0,6981 

- 20,7 

252,5 

619,19 

0,4994 

0,5699 

0,8762 

19.000 

5,79 

190,00 

485,48 

0,6751 

- 22,6 

250,5 

616,76 

0,4791 

0,5511 

0,8694 

20.000 

6,10 

200,00 

465,63 

0,6527 

- 24,6 

248,5 

614,32 

0,4595 

0,5328 

0,8625 

21.000 

6,40 

210,00 

446,45 

0,6308 

- 26,6 

246,5 

611,86 

0,4406 

0,5150 

0,8556 

22.000 

6,71 

220,00 

427,91 

0,6095 

- 28,6 

244,6 

609,40 

0,4223 

0,4976 

0,8487 

23.000 

7,01 

230,00 

410,01 

0,5888 

- 30,6 

242,6 

606,93 

0,4046 

0,4807 

0,8419 

24.000 

7,32 

240,00 

392,71 

0,5686 

- 32,5 

240,6 

604,44 

0,3876 

0,4642 

0,8350 

25.000 

7,62 

250,00 

376,01 

0,5489 

- 34,5 

238,6 

601,95 

0,3711 

0,4481 

0,8281 

26.000 

7,92 

260,00 

359,89 

0,5298 

- 36,5 

236,6 

599,45 

0,3552 

0,4325 

0,8212 

27.000 

8,23 

270,00 

344,33 

0,5112 

- 38,5 

234,7 

596,93 

0,3398 

0,4173 

0,8144 

28.000 

8,53 

280,00 

329,32 

0,4931 

- 40,5 

232,7 

594,41 

0,3250 

0,4025 

0,8075 

29.000 

8,84 

290,00 

314,85 

0,4754 

- 42,5 

230,7 

591,87 

0,3107 

0,3881 

0,8006 

30.000 

9,14 

300,00 

300,90 

0,4583 

- 44,4 

228,7 

589,32 

0,2970 

0,3741 

0,7937 

31.000 

9,45 

310,00 

287,45 

0,4417 

- 46,4 

226,7 

586,76 

0,2837 

0,3605 

0,7869 

32.000 

9,75 

320,00 

274,49 

0,4255 

- 48,4 

224,8 

584,20 

0,2709 

0,3473 

0,7800 

33.000 

10,06 

330,00 

262,01 

0,4097 

- 50,4 

222,8 

581,61 

0,2586 

0,3345 

0,7731 

34.000 

10,36 

340,00 

249,99 

0,3944 

- 52,4 

220,8 

579,02 

0,2467 

0,3220 

0,7662 

35.000 

10,67 

350,00 

238,42 

0,3796 

- 54,3 

218,8 

576,42 

0,2353 

0,3099 

0,7594 

36.000 

10,97 

360,00 

227,29 

0,3652 

- 56,3 

216,8 

573,80 

0,2243 

0,2981 

0,7525 

36089 

11 

361 

226,32 

0,3639 

- 56,5 

216,65 

573,57 

0,2234 

0,2971 

0,75187 

37.000 

11,28 

370,00 

216,627 

0,3483 

- 56,5 

216,6 

573,57 

0,2138 

0,2844 

0,7519 

38.000 

11,58 

380,00 

206,462 

0,3320 

- 56,5 

216,6 

573,57 

0,2038 

0,2710 

0,7519 

39.000 

11,89 

390,00 

196,773 

0,3164 

- 56,5 

216,6 

573,57 

0,1942 

0,2583 

0,7519 

40.000 

12,19 

400,00 

187,539 

0,3016 

- 56,5 

216,6 

573,57 

0,1851 

0,2462 

0,7519 

41.000 

12,50 

410,00 

178,739 

0,2874 

- 56,5 

216,6 

573,57 

0,1764 

0,2346 

0,7519 

42.000 

12,80 

420,00 

170,351 

0,2739 

- 56,5 

216,6 

573,57 

0,1681 

0,2236 

0,7519 

43.000 

13,11 

430,00 

162,357 

0,2611 

- 56,5 

216,6 

573,57 

0,1602 

0,2131 

0,7519 

44.000 

13,41 

440,00 

154,738 

0,2488 

- 56,5 

216,6 

573,57 

0,1527 

0,2031 

0,7519 

45.000 

13,72 

450,00 

147,477 

0,2371 

- 56,5 

216,6 

573,57 

0,1455 

0,1936 

0,7519 

46.000 

14,02 

460,00 

140,556 

0,2260 

- 56,5 

216,6 

573,57 

0,1387 

0,1845 

0,7519 

47.000 

14,33 

470,00 

133,960 

0,2154 

- 56,5 

216,6 

573,57 

0,1322 

0,1758 

0,7519 

48.000 

14,63 

480,00 

127,674 

0,2053 

- 56,5 

216,6 

573,57 

0,1260 

0,1676 

0,7519 

49.000 

14,94 

490,00 

121,683 

0,1957 

- 56,5 

216,6 

573,57 

0,1201 

0,1597 

0,7519 

50.000 

15,24 

500,00 

115,973 

0,1865 

- 56,5 

216,6 

573,57 

0,1145 

0,1522 

0,7519 

51.000 

15,54 

510,00 

110,530 

0,1777 

- 56,5 

216,65 

573,57 

0,1091 

0,1451 

0,75187 

52.000 

15,85 

520,00 

105,344 

0,1694 

- 56,5 

216,6 

573,57 

0,1040 

0,1383 

0,7519 

53.000 

16,15 

530,00 

100,400 

0,1614 

- 56,5 

216,6 

573,57 

0,0991 

0,1318 

0,7519 

54.000 

16,46 

540 . 0 Ü 

95,689 

0,1539 

- 56,5 

216,6 

573,57 

0,0944 

0,1256 

0,7519 

55.000 

16,76 

550,00 

91,198 

0,1466 

- 56,5 

216,6 

573,57 

0,0900 

0,1197 

0,7519 

56.000 

17,07 

560,00 

86,919 

0,1398 

- 56,5 

216,6 

573,57 

0,0858 

0,1141 

0,7519 

57.000 

17,37 

570,00 

82,840 

0,1332 

- 56,5 

216,6 

573,57 

0,0818 

0,1087 

0,7519 

58.000 

17,68 

580,00 

78,952 

0,1270 

- 56,5 

216,6 

573,57 

0,0779 

0,1036 

0,7519 

59.000 

17,98 

590,00 

75,247 

0,1210 

- 56,5 

216,6 

573,57 

0,0743 

0,0988 

0,7519 

60.000 

18,29 

600,00 

71,716 

0,1153 

- 56,5 

216,6 

573,57 

0,0708 

0,0941 

0,7519 

61.000 

18,59 

610,00 

68,351 

0,1099 

- 56,5 

216,6 

573,57 

0,0675 

0,0897 

0,7519 

62.000 

18,90 

620,00 

65,143 

0,1047 

- 56,5 

216,6 

573,57 

0,0643 

0,0855 

0,7519 

63.000 

19,20 

630,00 

62,087 

0,0998 

- 56,5 

216,6 

573,57 

0,0613 

0,0815 

0,7519 

64.000 

19,51 

640,00 

59,173 

0,0951 

- 56,5 

216,6 

573,57 

0,0584 

0,0777 

0,7519 

65.000 

19,81 

650,00 

56,396 

0,0907 

- 56,5 

216,6 

573,57 

0,0557 

0,0740 

0,7519 
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5.2 Ermittlung des relativen Fehlers in Prozent 


Mittels Verweises auf Tabelle 5.2 und Gl. (3.7) können nun die Machzahlen in Abhängigkeit 
der Höhe von oben nach unten zunehmend mit einem Intervall von 1000 ft und der kalibrier¬ 
ten Geschwindigkeit von li nk s nach rechts im Intervall von 10 kt wie in Tabelle 5.3 angege¬ 
ben ermittelt werden: 


Tabelle 5.3 Machzahlen in Abhängigkeit von H und CAS 



Aus den Werten der Tabelle 5.3 sowie Gl. (4.12) wird wiederum jeweils die „True Airtempe- 
rature“ berechnet (Tabelle 5.4). 


Tabelle 5.4 TAT in Kelvin in Abhängigkeit von H und CAS 
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Abgeleitet aus der Tabelle 5.4 und Gl. (4.19) können die gesuchten TAS-Werte berechnet 
werden (Tabelle 5.5). 


Tabelle 5.5 TAS in kt in Abhängigkeit von H und CAS 



Folgend wird mit Tabelle 5.6 die wahre Geschwindigkeit erneut über die zu prüfende Faust¬ 
formel aus Gl. (4.20) berechnet, wobei die Ergebnisse der Faustformel den Index ROT tragen: 


Tabelle 5.6 TASrot in kt über Faustformel in Abhängigkeit von H und CAS 
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Werden die Abweichungen der Ergebnisse aus Tabelle 5.5 und Tabelle 5.6 in einer weiteren 
Tabelle über Gl. (4.22) in Prozent dargestellt, so ergibt sich der relative Fehler in Abhängig¬ 
keit von Höhe und kalibrierter Geschwindigkeit in Prozent (Tabelle 5.7). 


Tabelle 5.7 Relativer Fehler der Faustformel in Abhängigkeit von H und CAS 



Um die Ergebnisse aus Tabelle 5.7 besser auswerten zu können, wird die gesamte Tabelle mit 
einer bedingten Formatierung versehen, wodurch die Abweichungen farblich akzentuiert wer¬ 
den. Hierbei werden Abweichungen, die kleiner als 5 % sind in Grün dargestellt. Größere 
Abweichungen nehmen von Gelb zu Rot zu (Tabelle 5.8). 
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5.3 Faustformel für geringe Geschwindigkeiten und Höhen 

Die Faustformel aus Gl. (1.6) darf bei hohen Geschwindigkeiten und Höhen unter Reiseflug¬ 
bedingungen angewendet werden. Der vollständigkeitshalber wird an dieser Stelle eine weite¬ 
re Faustformel mit Gl. (5.1) vorgestellt, die für niedrige Geschwindigkeiten und Höhen ver¬ 
wendet werden kann und somit für kleinmotorige Luftfahrzeuge angedacht ist (z.B. Cessna, 
Piper). 

Der Fahrtmesser berücksichtigt Dichte- bzw. Temperaturänderungen mit steigender Höhe 
nicht. Daher gilt bei steigender Höhe über Meereshöhe, das IAS kleiner als TAS wird (siehe 
Kapitel 2). Ferner gilt im Hochgeschwindigkeitsbereich (Reiseflug, Klappen eingefahren), 
dass IAS ungefähr CAS entspricht. Der Fahrmesser zeigt pro 1000 ft ca. 2 % zu wenig an. 

(Kühr 1987) 

Beispiel nach Kühr 1987: 

IAS: 100 kt 

Flughöhe: 5000 ft 
TAS: ? 

Lösung: 

2% -5 =10% (0,1) 

100 - 0,1 =10 
100+10 =110 
—> TAS =110 kt 

Die Faustfonnel nach Kühr 1987 lässt sich somit wie in Gl. (5.1) darstellen. 


v ROT n 


v c + 2 % ■ v c ■ 


H 

1000 /t 


(5.1) 


Vrot,, TAS aus Faustformel für geringe Geschwindigkeiten und Höhen 

Nun wird in Excel der relative Fehler der Faustformel in Prozent aus Gl. (5.1) nach gleichem 
Schema wie bereits bei der vorangegangenen Faustformel aus Gl. (1.6) ermittelt und an¬ 
schließend wieder mittels bedingter Formatierung farblich der Gültigkeitsbereich hervorgeho¬ 
ben, wobei Abweichungen kleiner gleich 5 % Grün dargestellt werden und die Abweichungen 
von Grün über Gelb nach Rot kontinuierlich zunehmen (Tabelle 5.9). 
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Tabelle 5.9 Relativer Fehler für Gl. (5.1) für geringe Höhen und Geschwindigkeiten 
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6 Matlab Live Skript 

6.1 3D-Diagramm 

Bislang wurde davon ausgegangen, dass A7=0, also folglich ISA-Bedingungen vorzufinden 
sind. Nun soll darüber hinaus geprüft werden, ob die Faustfonnel auch gültig bleibt, wenn die 
Referenztemperatur von den ISA-Bedingungen abweicht. Dies führt dazu, dass der relative 
Fehler von den 3 Parametern CAS, H und A T abhängig wird. Somit können die Ergebnisse 
nicht mehr 2-dimensional dargesteht werden. Daher wird über die vorangegangenen Glei¬ 
chungen ein Matlab-Live-Skript geschrieben, um die Abweichungen nun als relativen Fehler 
in Prozent mit vorgegebenen Anfangs- und Endwerten in festgelegten Intervallen zu errech¬ 
nen. 


6.2 Aufbau des Live Skripts 


Zuerst werden relevante Konstanten definiert (Bild 6.1). 


Constants: 

g 

= 

32.17404856; 

% [ft / s A 2] 

T0 

= 

288.15; 

% [K] 

H_tnopopause 

= 

36089.23885; 

% [ft] 

K_A 

= 

5.25587743244413; 


gamma 

= 

1.4; 


R_SI 

= 

287.053; 

% [(m A 2 / s A 2) / K] 

R 

= 

3089.81277682306; 

% [ (ft A 2 / s A 2) / K] 

L 

= 

0.0019812; 

% [K / ft] 

m_s_kt 

= 

1.9438444941; 

% factor to convert [m/s] -> [kt] 


Bild 6.1 Definierte Konstanten im Matlab-Code 


Um bei den Anfangs- und Endwerten sowie Intervallen flexibel zu bleiben, werden diese als 
Dropdown Menüs integriert (Bild 6.2). 
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ränge for variables 


dT_Start 

-10 

▼ ; dT_End 

15 


H_Start 

(o 

- ; H_End 

65000 

▼ 

JJ 

CAS_Start = 

10 

▼ ; CAS_End = 

400 

■w 

y 


Interval AT, H, CAS 


interval_dT= 

(0.2 

▼ 

-- 

interval_H= 

100 

▼ 

- j 

interval_CAS= 

Gl 

• 


Bild 6.2 Start- und Endwerte sowie Intervalle als Dropdownmenü 


Da zunächst noch nicht klar ist, wie die Ergebnisse graphisch aussehen werden, werden die 
Bezugswerte für die Farben und Punktgrößen auch via Dropdownmenüs eingepflegt, um diese 
gegebenenfalls anpassen zu können. Zudem werden Lösungsmengen definiert, auf die im spä¬ 
teren Code verwiesen wird, um im Diagramm später selektieren zu können, welcher Teil der 
Lösungsmenge angezeigt werden, respektive ausgeblendet soll und um die einzelnen Bereiche 
dann zur visuellen Unterscheidung die Farben zuzuweisen. Die finalen Lösungsmengen wur¬ 
den mit den berechneten TAS-Werten aus der Faustformel in Abhängigkeit der Parameter 
CAS, H und A T bezüglich der Abweichungen wie folgt definiert (Bild 6.3): 


Abweichungen 
bis 2 % 

2 % ... 5 % 

5 % ... 10 % 

10 % ... 15 % 
>15 % 


Dunkelgrün 

Grün 

Gelb 

Orange 

Rot. 
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graphic dotsizes 

dotsize0=0.1; 

dotsizel=[0.5 

- 

dotsize2=( 0.5 

- 

dotsize3=(o.2 

- 

dotsize4=( 0.5 

- 

dotsize5=[a2 

- 


dot colors and linked conditions 


cl= 

' darkgreen 

- 

11= 

z 


solution_setl= 

■gi 


■; % 

>= E 



c2= 

'[ green 

» 

12= 

GlI 

W > 

solution_set2= 

1v 


% 

>= £ 

> 

ii 

c3= 

yellow 

- 

13= 

Qo 

- ; 

solution_set3= 

'(y 


I'; % 

>= £ 

> 

12 

c4= 

'[ orange 

V J 

14= 


T y 

solution_set4= 

■|v 


l'j % 

>= £ 

> 

13 

c5= 

dark red 

w J 

15= 

14 ; 


solution_set5= 



% 

< £ 




Bild 6.3 Punktgrößen, Farben, Lösungsbereiche als Dropdownmenü 

Über den Matlab-Code (grauer Hintergrund) werden die verwendeten Formeln jeweils via 
Formeleditor dargestellt (weißer Hintergrund), um entsprechende Rechenschritte besser und 
schneller nachvollziehen zu können. Zudem werden rechts jeweils in entsprechender Zeile 
kurze Erklärungen (grüner Text) zu dem Abschnitt angegeben (Bild 6.4). 

Da es üblich ist den Code selbst und Definitionen in englischer Sprache zu formulieren, wird 
dies für das gesamte Matlab-Skript jeweils konsequent beibehalten, obgleich diese Projektar¬ 
beit selbst in deutscher Sprache verfasst wurde. 

Somit kann der Code bei Interesse/Bedarf im originalen Matlab-File leicht nachvollzogen 
werden und es werden im Rahmen dieser schriftlichen Ausarbeitung nur die entsprechenden 
Ergebnisse verarbeitet, ohne im Detail auf den Code einzugehen. 

Es werden insgesamt 2 Matlab-Files geschrieben. Im ersten File werden die Lösungsmengen 
rein farblich wie zuvor beschrieben unterschieden. Das Ergebnis ist somit ein dreidimensiona¬ 
les Diagramm in Abhängigkeit von den 3 Parametern CAS, H und A T. Im zweiten File wer¬ 
den die Abweichungen über die 3. Achse als relativer Fehler auch als Zahlenwert in Prozent 
entsprechend angezeigt. Die Lösungsmenge(n) ist daher eine Fläche in Abhängigkeit von den 
2 Parametern CAS und H. Die farblichen Akzente bleiben identisch. 
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Calculation of the relative error e for each possile combination of dT, CAS, H within the chosen ränge 

for k=dT_Start:interval_dT:dT_End %range for dT and Start of loop 1 
dT=k; %define dT for next loop 3 (inner loop to create the tab) 

for m=H_Start:interval_H:H_End %range for H and Start of loop 2 
H=m; %define H for loop 3 

for n=CAS_Start:interval_CAS:CAS_End %range for CAS and Start of loop 3 
CAS =n; %define CAS for each loop 3 (n/10 times per loop 1) 


Troposphäre: 

a - x r .R.T = v y • R • (T 0 + A T — I. • H) 

if H<= H_tropopause; a = sqrt( gamma * R_SI *(T0+dT-L*H))* m_s_kt ; % result troposphere [kt] 

Stratosphäre: 

a = v y •«•'/, = yj-f • R • T 0 + AT - l. • Hj) 

eiseif H > H_tropopause 

a = sqrt( gamma * R_SI * ( T0 + dT - L * H_tropopause ) ) * m_s_kt ; % result [kt] stratosphere 

end 


Troposphäre: 


R»L 


i=(\ —L»g.V 

V ?„ + A7-J 


if H<= H_tropopause; DELTA = (1 - ((L * H) / (T0 + dT))) A K_A; 

Stratosphäre: 


Tr = T 0 +AT 


K„ 


x 

R • T t 


6 = Sj • e 


~ h t ) 


&i = 



\ K * 

T„ + A T) 


eiseif H > H_tropopause 

DELTA_tropopause = (1 - (L * H_tropopause) / (T0 + dT)) A (K_A) ; 
T_tropopause = T0 + dT - L*H_tropopause; 

K_B = g/(T_tropopause * R); % g/(R*TT) -> [1/ft] 

DELTA = DELTA_tropopause * exp((-K_B)*(H-H_tropopause)); 

end 


Bild 6.4 


Auszug aus dem Matlab-Code 
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7 Auswertung der Ergebnisse 

7.1 ISA-Tabellen 

Zunächst wird die Ergebnisstabelle (Tabelle 7.1) weiterverarbeitet. Von li nk s nach rechts mit 
steigenden CAS-Werten und von oben nach unten mit steigender Höhe //wird jeweils die re¬ 
lative Abweichung in Prozent angegeben, die sich aus der Faustformel zu den Ergebnissen 
aus der Berechnung mit den Gleichungen aus der genutzten Literatur aus Kapitel 2 bis 4 erge¬ 
ben. Abweichungen bis 5 % sind in Grün zu sehen, während größere Abweichungen von Gelb 
zu Rot kontinuierlich zunehmen. 



Da die zu prüfende Faustfonnel aus Gl. (1.6) für Luftfahrzeuge unter Reisebedingungen mit 
hohen Geschwindigkeiten und Höhen angedacht ist, wird in einer weiteren Tabelle der Reise¬ 
flugbereich von gängigen zivilen Flugzeugen in Abhängigkeit von der Machzahl dargestellt 
(Tabelle 7.2). Hierbei wird davon ausgegangen, dass je nach Antriebsleistung und Variante 
(Propeller, Jet, Turboprob, Turbofan, etc.) die Reisfluggeschwindigkeit zwischen 0,70 Mach 
und 0,85 Mach liegt. Der entsprechende Bereich ist Schwarz hinterlegt. Größere Machzahlen 
als 0,85 Mach sind mit rotem Hintergrund und kleinere Machzahlen als 0,7 Mach mit gelbem 
Hintergrund dargestellt. Wie bereits in Tabelle 7.1 werden die Machzahlen analog in Abhän¬ 
gigkeit von den Parametern CAS und H angegeben. 
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Tabelle 7.2 Reiseflugbedingungen (schwarz) 

H IcAS [kt] 
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Legt man nun beide vorangegangenen Tabellen übereinander (Tabelle 7.3), so ist deutlich zu 
sehen, dass der Machzahlenbereich unter Reiseflugbedingungen (schwarz-transparent) genau 
in den Bereich fällt, bei dem die Abweichungen der Faustformel unter 5 % hegen (Grün). 

Tabelle 7.3 Überlagerung „cruise speed“ und relativer Fehler 


CAS 



Wird nun der Reiseflugbereich eines zivilen Flugzeugs weiter eingegrenzt und in Bezug auf 
die Abweichungen der Faustformel hin betrachtet (Tabelle 7.4), so ist noch deutlicher zu er¬ 
kennen, dass die Faustformel im gesamten Bereich unter Reisflugbedingungen bedenkenlos 
angewendet werden darf, ohne dass Abweichungen größer als 5 % resultieren. 
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Tabelle 7.4 Anwendungsbereich der Faustformel aus Gl. (1.6) 
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7.2 3D-Diagramm in Matlab 


Das Matlab-Skript ergibt im Abhängikeit der Parameters A T, CAS und H nun ein dreidimen¬ 
sionales Diagramm, dass wie in Kapitel 6.2 beschrieben aufgebaut ist. Wenngleich das Re¬ 
sultat annähernd einem Volumen gleicht, so werden eine Vielzahl einzelner Ergebnis-Punkte 
dargestellt, bei denen die Intervalle und Anfangs- und Endwerte so lange angepasst wurden, 
bis sich die visuell ansprechendsten Ergebnisse ergaben. 

Die Abgrenzungen, die auch in der Legende jeweils angeben werden, sind an dieser Stelle 
noch einmal aufgeführt: 


Abweichungen: 
bis 2 % 

2 % ... 5 % 

5 % ... 10 % 

10 % ... 15 % 
>15 % 


Dunkelgrün 

Grün 

Gelb 

Orange 

Rot. 


Nachfolgend werden mit Bild 7.1 bis Bild 7.5 das Ergebnis aus verschieden Perspektive dar¬ 
gestellt. 



Bild 7.1 Abweichung der Faustformel aus Gl. (1.6), Perspektive 1 
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Bild 7.2 Abweichung der Faustformel aus Gl. (1.6), Perspektive 2 



Bild 7.3 Abweichung der Faustformel aus Gl. (1.6), Perspektive 3 
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CAS [kt] 400 -10 dT[°C] 



E < 

• 

2 

• 

5 

• 

10 

• 

15 


£ > 

• 

15 


Bild 7.4 Abweichung der Faustformel aus Gl. (1.6), Perspektive 4 


x 10 4 




E < 

• 

2 

• 

5 

• 

10 

• 

15 


E > 

• 

15 


Bild 7.5 


Abweichung der Faustformel aus Gl. (1.6), Perspektive 5 
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Da schlussendlich die Lösungsmenge, bei der die Abweichungen kleiner als 5 %, bzw. kleiner 
2 % von besonderem Interesse sind, werden in den folgenden Darstellungen von Bild 7.7 bis 
Bild 7.11 alle Werte mit Abweichungen größer als 5 % über das eingangs angelegte 
Dropdownmenü (Bild 7.6) ausgeblendet. 


dot colors and linked conditions 


cl= 

dark green 

v'j 

11 = 

2 

«■ ; 

solution_setl= 

■fv 

- 'j % 

>= £ 



c2= 

green 

- 

12= 

> 

- ; 

solution set2= 

' [y 

- % 

>= £ 

> 

11 

c3= 

' yellow 

* 

13= 

|10 

- ; 

solution_set3= 

i n 

- % 

>= £ 

> 

12 

c4= 

orange 

» 

14= 

J 5 

- ; 

solution_set4= 

i n 

’ 'i % 

>= £ 

> 

13 

c5= 

dark red 

- 

15= 

14 ; 


solution set5= 

n 

zhJL 

< £ 




Bild 7.6 Dropdownmenü für Lösungsmengen und Farben 



Bild 7.7 


Abweichung kleiner 5 %, Perspektive 1 






































H [ft] x H [ft] 
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x 10 

0 , 



100 


200 


300 


CAS [kt] 

400 

- 10 dT [°C] 

7.8 Abweichung kleiner 5 %, Perspektive 2 


E < 

• 2 
• 5 


10 4 

0 -i 



Bild 7.9 


Abweichung kleiner 5 %, Perspektive 3 








56 



E < 

• 2 
• 5 


15 


Bild 7.10 Abweichung kleiner 5 %, Perspektive 4 
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Bild 7.11 


Abweichung kleiner 5 %, Perspektive 5 
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Mit Blick auf die Tabelle 7.2 und den Abgleich mit den Ergebnissen aus Matlab, ist zu er¬ 
kennen, dass A T zwar bei negativen A T- Werten den Gültigkeitsbereich in niedrigere Höhen, 
sowie bei positiven Ar-Werten in höhere Höhen verschiebt, im geprüften Bereich zwischen - 
10 °C und +15 °C liegen die sich ergebenen TAS-Werte der Faustfonnel allerdings unter 
Reiseflugbedingungen noch bei Abweichungen unter 5 %. Dies ist exemplarisch für die AT- 
Werte bei - 10 °C und + 15 °C nach gleichem Schema wie in Tabelle 7.3 jeweils in Tabelle 
7.5 und Tabelle 7.6 aus Excel aufgezeigt. Zu sehen sind in Schwarz-transparent der Reise¬ 
flugbereich und in Grün Abweichungen der Faustformel unter 5 %. 


Tabelle 7.5 Anwendungsbereich der Faustformel aus Gl. (1.6) bei A7=-10 °C 



Tabelle 7.6 Anwendungsbereich der Faustformel aus Gl. (1.6) bei A7=15 °C 



Relevanter Bereich (Reiseflughöhe) 
H = 30000 - 45000 ft 
CAS = 200 - 320 kt 
M= 0,7-0,85 


CAS [kt] 


Bei Temperaturabweichung, die noch weiter von den Standardbedingungen nach ISA abwei- 
chen, verschiebt sich der Bereich entsprechend weiter nach unten bei AT-Werten größer als 
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15 °C, respektive nach oben bei AT-Werten kleiner als -10 °C. Dabei ergeben sich rechne¬ 
risch dann Abweichungen zwischen 5 % und 10 %, was aber immer noch eine recht gute An¬ 
näherung darstellt, zumal die Faustfonnel sehr schnell und ohne großen Aufwand entspre¬ 
chende Ergebnisse liefert. 


7.3 3D-Diagramm mit \7=0 

Wird nun über das 3D-Diagramm aus Matlab der AT-Wert und das Intervall auf null gesetzt, 
ergibt sich wieder ein zweidimensionales Diagramm in Abhängigkeit der 2 Parameter CAS 
und H unter ISA-Bedingungen (Bild 7.12) und kann somit direkt mit der Tabelle 7.1 aus 
Excel verglichen werden. 



Zum Abgleich wird das Ergebnis aus Bild 7.12 in der Seitenansicht dargestellt (Bild 7.13). 
Werden nun beide Ergebnisse unmittelbar übereinander dargestellt (Bild 7.14), so sieht man 
wie zu erwarten, dass die Ergebnisse in Matlab die gleichen Resultate liefert, wie zuvor be¬ 
reits die Tabellen aus Excel über bedingte Formatierungen. 
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Bild 7.13 3D-Diagramm mit A7 und Intervall von 0, Seitenansicht 
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Bild 7.14 Vergleich ISA-Tabelle und Matlab-Diagramm 
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7.4 2D-Diagramm unter ISA-Bedingungen in Matlab 

Eine weitere Möglichkeit die Ergebnisse in Matlab unter ISA-Bedingungen, d.h. mit einem 
A T von null darzustellen, bietet eine Darstellung über die beiden variablen Parameter CAS 
und H auf jeweils einer Achse und dem Ergebnis als relativer Fehler in Prozent auf der dritten 
Achse. Bei kleineren Abfrageintervallen sind noch einzelne Ergebnispunkte zu erkennen 
(Bild 7.15) Das Ergebnis liefert annähernd eine Fläche, sofern die Intervalle eng genug ange¬ 
setzt werden. Zu bedenken ist, dass sehr große Matrizen verarbeitet werden, die entsprechen¬ 
den Rechenaufwand fordern und nur unter Rechnern mit hohem Leistungsvermögen in an¬ 
nehmbarer Zeit Ergebnisse ausgegeben werden. 



In Bild 7.16 bis Bild 7.18 wird das Ergebnis aus dem Matlab-Skript jeweils in verschiedenen 
Perspektiven mit kleineren Intervallen dargesteht. 

In der Draufsicht (Bild 7.19) ergibt das Diagramm das bekannte Ergebnis aus Bild 7.14. 
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Bild 7.16 2D-Diagramm, Perspektive 1 



Bild 7.17 


2D-Diagramm, Perspektive 2 
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40% 



CAS [kt] 


Bild 7.18 2D-Diagramm, Perspektive 3 



Bild 7.19 


2D-Diagramm, Draufsicht 
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8 Zusammenfassung 

Es wurden in dieser Projektarbeit neben allgemeinen Zusammenhängen aus dem Fachgebiet 
der Flugmechanik das „Comressibility Correction Chart,“ diverse Gleichungen aus verschie¬ 
denen Fiteraturquellen sowie eine, bzw. zwei Faustformel vorgestellt. Dabei können die vor¬ 
gestellten Methoden allesamt zur Umrechnung von CAS in TAS genutzt werden, wobei die 
Faustformel mittels Excel und Matlab geprüft wurde. Es ist im Rahmen der Projektarbeit ge¬ 
lungen einen mathematischen Zusammenhang zwischen der zu prüfenden Faustformel aus Gl. 
(1.6) und den bekannten Gleichungen der Flugmechanik herzustellen, so dass daraus eine re¬ 
lative Abweichung in Prozent über definierte Geschwindigkeits- und Höhenbereiche ennittelt 
werden konnte. Die relativen Fehler wurden in Exceltabellen in festgelegten Intervallen er¬ 
rechnet und anschließend über bedingte Formatierungen farblich fließend von Grün nach Rot 
der Abweichung entsprechend visualisiert. Dabei zeigt sich deutlich, dass die besagte Faust¬ 
formel hinreichend genaue Ergebnisse mit Abweichung unter 5 % liefert und somit ohne Be¬ 
denken, z.B. im Blindflug unter Reiseflugbedingungen genutzt werden kann. Zudem wurde in 
Matlab mit der Temperatur ein weiterer Parameter eingeführt und es hat sich gezeigt, dass die 
Faustformel auch bei Abweichungen von ±15° unter Reisflugbedingungen hinreichend ge¬ 
naue Ergebnisse erbringt. 
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9 Ausblick 

Wenngleich im Rahmen dieser Projektarbeit lediglich eine, respektive zwei Faustformeln ge¬ 
prüft wurden, so zeigen die Ergebnisse, welch gute Annährung eine Faustformel auf die exak¬ 
ten Ergebnisse liefern kann. Zudem werden Mitteln und Wege übermittelt, wie Faustformeln 
mit einfachen Werkzeugen, wie Excel mit relativ bescheidenem Aufwand über größere Berei¬ 
che geprüft werden können. Daher ist diese Projektarbeit nicht ausschließlich als Verifikation 
einer einzelnen Faustformel zu verstehen, sondern vielmehr als exemplarischer Nachweis da¬ 
für, dass aus einer Faustformel recht genaue Ergebnisse resultieren können, ohne dafür tech¬ 
nische Hilfsmittel nutzen zu müssen. Darüber hinaus kann das angewandte Schema auf weite¬ 
re Faustformeln übertragen werden, wodurch „neue“ Faustformeln leicht geprüft werden kön¬ 
nen, bevor man sie auf das eigene „Kniebrett“ überträgt. 
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Anhang Aufbau der ergänzenden Exceldatei 

Allgemeiner Aufbau: 

Die Excel-Datei besteht insgesamt aus 11 einzelnen Untertabellen, die folgend von links nach 
rechts dem Inhalt entsprechend kurz beschrieben werden: 

• „Konstanten“ 

In der Tabelle „Konstanten“ sind häufig genutzte Konstanten und Faktoren 
zum Konvertieren von Einheiten enthalten. 

• „ISA TAB blank“ 

In der Tabelle „ISA TAB bla nk “ wird zunächst eine ISA-Tabelle über die zu¬ 
vor beschrieben Gleichungen aus Kapitel 2 bis Kapitel 4 bis zu einer Höhe von 
65000 ft generiert und daraus abgeleitet die Machzahlen, TAT- und TAS- 
Werte bis zu Geschwindigkeiten von 400 kt errechnet (aus CAS). 

• „cruise speed DeltaT = -10 °C“ 

Gleiches Prozedere wie Tabelle „relative error colored“ mit einem AT-Wert 
von -10°C. 

• „cruise speed Delta T = +15 °C“ 

Gleiches Prozedere wie Tabelle „relative error colored“ mit einem AT-Wert 
von +15°C. 

• „cruise speed Delta T = +30 °C“ 

Gleiches Prozedere wie Tabelle „relative error colored“ mit einem AT-Wert 
von +30 °C. 

• “compressibility correction” 

Es wird ein CAS-Wert und eine Höhe H gewählt (Input, Grau hinterlegt), 
wodurch in dem enthaltenen Chart u.a. der AE-Wert angezeigt wird. Weiter 
rechts vom Chart sind die zu Grunde liegenden Tabellen, auf die im interakti¬ 
ven Chart verwiesen wird zu finden. Über den Tabellen sind die verwendeten 
Gleichungen abgebildet, auf denen die Tabellen basieren. 



“relative error” 

Neben den Tabellen aus „ISA TAB blank“, die übertragen wurden, werden die 
TAS-Werte erneut über die zu prüfende Faustformel berechnet und über die 
Abweichungen des relativen Fehlers in Prozent angeben. 

„relative error colored“ 

Die Tabelle „relativer error colored“ ist identisch mit der Tabelle „relative er¬ 
ror.“ Es wurden zusätzlich die prozentualen Abweichungen der Faustformel 
über eine bedingte Formatierung farblich akzentuiert. Von Grün (kleiner 5 % 
Abweichungen) nehmen die Abweichungen über Gelb zu Rot kontinuierlich 
zu. 


„Faustformel_2 (nach Kühr)“ 

„Faustforme 1_2 (nach Kühr)“ ist nach dem gleichen Schema wie „relative error 
colored“ aufgebaut. Es wurde hierbei die Faustformel für niedrige Geschwin¬ 
digkeiten und Höhen nach (Kühr 1987) geprüft und die Abweichungen ent¬ 
sprechend farblich visualisiert. 

„cruise speed“ 

„cruise speed“ ist ein Duplikat von „ISA TAB blank“ in der die Tabelle mit 
den enthaltenden Machzahlen in Abhängigkeit von Höhe und CAS-Werten 
farblich überarbeitet wurden, um den Bereich abzugrenzen, in dem Reisflug¬ 
bedingungen vorherreschen. 

„Überlagerung cruise speed“ 

In „Überlagerung cruise speed“ werden die beiden farblich überarbeiteten Er¬ 
gebnis-Tabellen „cruise speed“ und „relative error colored“ mittels Transpa¬ 
renzeffekts übereinandergelegt, um zu demonstrieren, dass die Ergebnisse der 
Faustformel im vorgesehen Bereich (Reiseflug) die genausten Ergebnisse lie¬ 
fern. 



